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ABSTRACT 


Uhls  report  presents  rerulte  on  the  synthesis  of  time -optimal  control 
ftor  a  second-order  nonlinear  system.  The  problem  is  to  determine  the 
feedback  control  as  a  function  of  the  state  variables.  This  is  equivalent 
to  finding  the  switching  locus,  since  the  tine -optimal  control  is  a  relay 
control.  The  nonlinear  system  is  represented  by  a  soft  spring  charactered 
by  the  Duff?  jg  equation  with  negative  nonlinear  term. 

The  d  main  of  controllability  is  described.  Ftor  the  case  that  the 
coefficient  of  the  nonlinear  term  is  sufficiently  small  a  procedure 
is  developed  vhlch  determines  the  maximum  number  of  switchings.  Also 
the  influence  of  damping  Is  considered. 

An  analytic  expression  of  the  feedback  contro1  >  derived  for  the 
conservative  case.  This  expression  contains  a  hyper-ellipt.  4  n  ral  A 
the  first  kind.  Since  this  integral  cannot  be  solved  in  j.osed  form 
estimate  ;  are  derived  to  show  the  existence  or  non-exi'  „ence  o  - 
switching  vjint  in  certain  regions  of  the  domain  of  .ontrolle  >lllty. 


CHAPTER  I 


INTRODUCTION 


This  thesis  nakes  extensive  reference  to  work  (lone  by  Markus  end  Lee 
CL  3.  lbs  principles  end  concepts  which  they  developed  needed  In  the 
following  chapters  will  be  presented  in  this  chapter. 

The  systen  under  consideration  is 

1  +  F(xyi)  -  n 


x  «  y 

f  -  -  P(**t)  +  *  (1) 

with  F(0, 0)  •  0t  T(x,j)  in  the  class  in  the  real  maber  plane  R 2,  u 
is  la  the  oonpact  interval  Ot  -1  <  u  <  1. 

Neat  of  the  investigation  in  Chapters  II  throng  IV  considers  a 
'jeneervatlvQ  systen  for  which  ■  F_(xfy)  5  0.  Then  in  equation  (1) 

"  J 

l*(x,y)  is  replaced  by  f(x).  But  any  physical  systen  has  daaplng, 

Fy(x,y)  £  0,  and  in  section  3»2  a  certain  Influence  of  the  daaplng  is 
discussed.  Therefore  in  this  chapter  the  theoreas  and  definitions  are 
stated  for  systens  of  the  form  (1)  with  nonsero  F  (x,y). 

w 

1.1  Controllability 

An  lapertant  concept  is  controllability  which  is  defined  for  nonlinear 
systens  in  [1,  Chapter  6.1].  This  definition  is  stated  here  for  second- 
order  systens. 


' 


-2- 


Defint*oa:  The  domin  G  of  null  controllability  for  the  control  process 

2 

(1)  in  P.  end  u  iu  D  in  defined  at  the  set  of  all  initial  points 
(xo,yo)  CR  ,  each  of  vhich  can  be  steered  to  the  origin  by  some  bounded 
measurable  controller  n(t)  €0  in  some  finite  time  duration.  If  G  contains 
an  open  neighborhood  of  the  origin,  then  the  system  la  said  to  be  locally 
controllable  near  the  origin. 

G  1b  connected  because  each  point  in  G  is  joined  to  the  origin  by  a 

2 

continuous  solution  curve  lying  entirely  in  G.  G  is  open  in  R  if  and 
only  if  it  contains  a  neighborhood  of  the  origin,  because  the  solutions 
of  this  differential  equation  depend  continuously  on  the  initial  conditions. 

In  [1,  Chapter  6.1]  a  theorem  is  proved  which  states  a  condition  under 
which  a  system  is  locally  controllable.  This  condition  for  the  second* 
order  system  under  consideration  is  that  a  certain  matrix  has  rank  2. 

Vor  system  (l)  this  matrix  la 
’  0  1 

1  F  (0,0) 

»  w 

Thun  for  this  systen  Independent  of  whether  damping  is  present  or  not 
the  domain  G  of  controllability  is  an  open  connected  subset  of  R^ 
containing  the  origin. 

1.2  Existence  of  Optimal  Control 

Another  theorem  needed  for  further  discussion  is  the  theorem  on  the 
existence  of  an  optLaal  control  [1,  Chapter  4.2].  Specified  for  second- 
order  autonomous  systems  it  reads: 

Theo^cr  1.2  Consider  the  nonlinear  proceas  in  R 


x  -  y 

y  *  -  F(x,y)  +  u 


(1) 


-3- 


f(x,y)  is  C1  in  R2,  a  €  0,  with  the  following  condition*: 

(1)  The  *et*  of  ell  initial  point*  (*o»P0)  sad  *11  target 
point*  (x^y^)  are  noneaqpty  and  conpact. 

(2)  The  control  restraint  integral  0  1*  nonempty  and 


(3)  The  state  constraint*  h*(x,y)  >  0,  i-1,..*,  rf  if 

2 

they  exist,  belong  to  the  class  C  in  R  . 

(k)  The  fully  J  of  adadssible  controllers  consists  of 

all  nsasttrable  fanctions  u(t)  on  variotis  tin*  interval* 
^,^3  c  such  that  each  n(t)  ha*  a  response 

(x(t),y(t))  on  Ct0*tx3  ateering  (x(t0),y(t0)M*0*y0> 
to  (x(t^),y(t^))  -  (x^)  end  n(t)  en» 


r(*(t),y(t)) 


(3)  The  cost  for  each  u  «  J  is 

S 

C(u)  -  *(x1»y1)  +  J  F°(x(t),j(t),u(t))  dt  ♦  as: 

to  *0^1 

shore  f°  is  C1  in  H3,  g(x,y)  and  y(x,y)  are  continuous  in  R2. 

If 

(a)  J  is  nonssfty. 

(b)  There  exists  a  uniforn  bound  x2(t)  +  y2(t)  <  N  on  [t0,t^3 
for  all  responses  (x(t),y(t))  to  u  €  J 

(c)  The  set  {^(xjy.u),  x,  y  |  «  €  fl}  i*  convex  in  R3  tar  each  (x,y), 
than  there  exists  an  opt  Inal  controller  m(t)  on  [t0,t^3  in  P  nininiting 
C(»). 

For  the  tlas-optlnal  problen  another  class  of  controllers  is  defined. 
Definition:  A  is  the  class  of  all  nsasurahle  controllers  u(t)  €  0  on 


finite  tins  intervals  [0,^3  ««eh  that  the  response  (x(t),  y(t))  of 


TEXT  NOT  REPRODUCIBLE 


/ 


(1)  defined  on  re,t  I  with  -n.>  «:1  <r-  pci’"*.  (r(O).  y(O) )  first 

reaches  th*  c-r i>r5 *A  t  t  . 


1.3  The  Octicjttl  '■3.r..rol 


The  ffcai  '  f  Mr*  vj  -  u*  ‘  v* 


:v?  *.  *  .jiz-  the  tine- 


ootlaai  feedback  cent  rt  jlu  r  'oj  r  ^-.e  ' '  ..  ".r '  .conditions  on 

?(x,y). 

A  controller  u<'t)  on  r  Hr  t*  w.i»o  *  t»  .  '*  in  ’  Is  called 

tine -optimal  if  for  er.  uf  (t)  in  f  or  '  2  .«  *  J.i*  that  t,  *  t,' .  ly 

tiie  naoiaum  principle  [1,  '’hspter  IthH  optical  -or;-*; roller  is  a 
maximal  controller.  That  mans  the  following  condition?  must  hold: 

(a)  There  exists  an  absolutely  ocnvjiuous.  novhore -vanishing 
adjoint  vector  (rjj(t),  t]?(t))  or.  '‘Cjt,  “}  such  tnat 
(x(t),  y(t)),  (r^t),  and  u(0  .vrtisf/  the 

Mami.lt  cnLar  ryster? 

«• 

.  —  ft  -  -  ~ 


y  "  ?zr 

•  2 

(b)  The  'foniltcninn 


H(Th.*Vx,r>u'>  *  T'i  /-r;:  *'*»»)  ”■  u 

is  maximized  with  respect  tn  u  ror  nlaor-t  all  t  on  [0,t^  j. 

(c)  The  ’’rnS ltoni»n  ic  greater  cr  cijurI  to  zero  and 

constant  for  o.ll  t  on  T f- , t ^  k lor..  *•.■>  optimal  frajec*tory. 

-1  if  TJ-(t  >  <  0 

Ia 

u(t)  -  ngn  nLXt)  “  0  if  r^{t)  •  C 

a  ir  n,(t)  >  o. 

V 


*  / 
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That  la,  tha  time-cptiaal  controller  la  a  relay  controller. 

Definition:  A  controller  a(t)  on  [0,t^3  In  A  la  called  a  relay  controller 
if  there  exlata  a  finite  maber  of  switching  tinea 

°",o<Tl<—<tk"tl 

auch  that  on  each  open  interval  i=l,...,k,  u(t)  la  constant 

and  equals  +1  or  -1,  and  u(t)  switches  values  on  successive  Intervals. 

Writing  the  Hamiltonian  systea  for  systen  (1)  one  obtains  with  the 
mrlnal  controller 

*  -  y 

f  »  -  F(x,y)  +  1  for  u  =  +  1  (1+) 

or 

*  -  y 

y  -  -  F(x,y)  -  1  for  u  =  -  1  (l_) 

and 

-  Tk  fx(x,y) 

Vz  “  *  Tl  +  ^2  VX,y)  ^ 

which  is  the  adjoint  systen.  Mote  that  the  Hanlltonlan  systea,  (1+)  and 
(2)  or  (1_)  and  (2),  is  homogeneous.  Since  (r^Ct),  r)2(t))  is  a  nonvanishing 
vector  solution  in  C1  on  [0,^]  of  (2)  ?)2( t)  has  only  a  finite  number  of 
zeros  which  arc  simple. 

By  this  maximum  principle  the  control  u(t)  is  only  determined  as 
a  function  of  tine.  The  feedback  controller  is  a  function  of  x  and  y 
which  rolves  the  synthesis  problem. 

1.4  The  Switching  Locus 

The  following  theorem  referred  to  as  the  "Theorem  of  Interlacing 
Zeros"  is  very  basic  in  establishing  the  properties  of  the  switching 
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locus  for  the  construction  of  the  feedback  controller  [1,  Chapter  7  *13* 
Theorem  1.4:  Let  u(t)  for  t  €  [0,t^3  in  A  be  a  maximal  controller  for 
system  (1),  and  let  (x(t),  y(t))  and  (f^(t),  rfg(t))  *he  corresponding 
response  c C  the  Hamiltonian  system,  (1)  and  (2),  on  [0,t^].  Let  t2 
be  time  cp  0  <  <  t2  <  t^,  then  four  assertions  hold: 


(1) 

If  V2^j)  ~ 

t?2(t2)  ■ 

'  0  and  if  y(x^) 

=  0,  then  y(r2)  *  0. 

(2) 

If 

=  VT25 

=  0  and  if  y(x^ 

)  /  0  then  y(r2)  /  0, 

but  there 

is  a  sero  of  y(t)  on  (x^ 

»t2)* 

(3) 

if  y(Tx)  *= 

JT(t2)  » 

o,  y(t)  ^  0  on 

(^.Tg),  and  if  V2(ri)  =  0 

then  tjgtTg 

)  *  o. 

(4) 

if  y(Tx)  “ 

y(r2)  - 

0,  y(t)  /  0  on 

(r1»Tg),  and  if  TJ2(tx)  /  0 

then  Tj^Tg 

)  /  0,  but  there  is  a  sero  of  rj2(t)  on  (TpXg). 

In  other  words:  Provided  the  seros  of  y(t)  are  isolated,  they  either  coin¬ 
cide  with  the  teros  of  r^ft),  or  no  sero  of  y(t)  is  a  aero  of  T)g(t),  but 
these  two  sets  of  seros  are  interlaced. 

Proof:  Assume  tj2(x1)  =  T72(t2)  =  0 

The  Hamiltoniar  is  y  -  »J2  F(x, y)  +  Irjgl  3  coiut.  >0  on  [0,^3 
Therefore  y^)  *  f^Cfg)  y(T2^  >  °«  3inc«  the  *#ro*  of  1^(0  *** 

simple  and  by  equation  (2)  tj^(t2)  <  °*  Thua  either  y(T^)  *  0 

if  and  only  if  y(r2)  *  0,  or  if  y^)  t  0,  then  y^)  y(r2)  <  0.  So  there 
is  a  sere  of  y(t)  on  (x^Tg).  This  proves  assertions  (1)  and  (2). 

Now  assume  y^)  “  y(T2)  =  0  and  y(t)  /  0  on  (x^Xg).  The  condition 
that  the  Hamiltonian  is  constant  on  [0,t^  3  yields 

y(Tx)  rj^)  «  Kt2)r»2(T2) 

With  y(t)  /  0  on  (t^,t2)  the  uniqueness  property  of  differential  equations 
requires  yZ(t)  +  £2(t)  /  0  on  Therefore  if  tj^x^  *  0  then 
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tlgCtg)  -  0.  If  y^)  -  y(r2)  *  0,  y(t)  /  0  on  (r^Tg)  and  rjg^)  j  0  then 
f^(t2)  /  0  hy  the  above  equation.  If  tj2(t)  vanlehea  nowhere  on 
then  y(r^)  y(t2)  >  0  which  is  impossible  since  and  t2  are  consecutive 
taros  of  y(t).  IMs  ccmpletec  the  proof. 

Corollary  l.frt  It  is  easy  to  show  that  at  a  switching  point:  t)2(t)  *  0 
if  y(r)  >  0  than  f^r)  <  0 

and  if  y(r)  <  0  then  ^2(t)  >  0 

At  t  the  second  equation  of  the  adjoint  system  (2)  is 
^(t)  «  -  ^(t) 

and  the  Hamiltonian  is 

f^(r)  l(r)  >  0 
or  fl2(r)  j(x)  <  0 

therefore  u(t)  switches  from  +1  to  -1  in  y  >  0  and  from  *1  to  +1  in  y  <  0, 
going  forward  in  time  along  an  optimal  trajectory. 

Markus  and  Less  describe  the  domain  of  controllability  and  prove  the 
existence  of  an  optimal  control  for  the  cases 

(a)  attractive  force,  x  F(x,0)  >  0  for  all  x  /  0 

(b)  repulsive  force  Fx(x,0)  <  0  for  all  x  /  0 

2 

with  moonegative  damping  7y(x,y)  >  0  in  R  .  A  similar  theorem  with  proof 
for  a  different  condition  on  F(x,y)  will  be  presented  in  Chapter  II. 
the  synthesis  problem  is  solved  hy  determining  the  switching  locus.  For 
e  secomd  order  system  it  may  be  described  as  a  function  of  x:  W(x). 
Definition:  The  twitching  locus  V  consists  of  the  points  (x,y)  at  which 
the  ■***—*  responses  switch  from  the  solution  family  of  (l+)  to  (l_)  or 


vice  versa,  the  origin  is  included  in  W 


[ 
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Therofore  the  trajectory  through  the  origin  (or  a  segment  of  it)  auat 
belong  to  W.  This  indicates  &  way  for  the  construction  of  W.  Due  to  the 
theorem  of  interlacing  zeros  and  its  corollary  a  part  of  the  solution  of  (1+) 
through  the  origin  in  the  fourth  quadrant  is  an  arc  of  the  switching  locus, 
W^.  The  solution  of  (1_)  through  the  origin  in  the  second  quadrant  as  a  part 
of  the  switching  locus  is  called  w\  Starting  with  an  initial  point  on 

(or  *£)  and  integrating  backwards  in  time  until  r»2(^)  =  0  results  in  a 

r  2 

point  of  the  arc  v£  (or  W+)  of  the  switching  locus.  The  initial  values  on 

V*  (or  are  T72(0)  »  0,  ^(0)  =  -1  (or  tj^O)  =  +1),  the  sign  is  given 

by  corollary  1.4,  the  magnitude  can  be  chosen  at  will  since  the  system 

112 

is  homogeneous,  x(0),  y(0)  are  given  as  coordinates  of  W+  (or  W  ).  W 

2  O  *3  —  — 

(and  W+)  are  reflected  the  same  way  to  obtain  vr  (and  Vr),  and  so  on.  Ihe 
switching  locus  is  then  the  union  of  all  sets  W*  and  w\  where  k  can  be 
finite  or  infinite.  In  general  the  integration  is  performed  numerically 
so  that  the  switching  locus  is  determined  pointvise. 

In  the  case  of  attractive  force  Markus  and  Lee  have  shown  that  k  >  1, 
if  there  is  no  friction,  F  5  0,  k  =  ».  In  the  case  of  repulsive  force 
k  =  1  and  the  switching  locus  coincides  with  the  trajectory  through  the 
origin  in  the  second  and  fourth  quadrant  over  its  entire  length. 

In  parts  of  their  discussion  they  assume  that  w(x)  is  single-valued 
which  need  not  be  true. 

As  an  example  of  a  conservative  system  with  attractive  force  they 

construct  the  switching  locus  for  the  time-optimal  control  of  a  hard 

spring  (Duffing  equation  ).  As  can  be  seen  from  the  results  the  witching 

locus  W(x)  *  U  (V^  U  W*)  is  not  single-valued  for  k  >  3. 
fc=l 

They  also  present  some  Interesting  results-mainly  computer  solutions- 


\ 


for  the  notion  of  a  forced  pendulum 


3— ark:  Later  in  the  analysis  singular  points  are  discussed.  These  are 
points  in  the  (x,y)  -  plane  for  which  x  and  $  are  zero  simultaneously. 

For  system  (1)  a  singular  point  can  exist  only  on  the  x-axis  with  x- 
coordlnate  a  real  root  of  F(x,0)  -  u  =  0  for  u  =  +1  or  -1.  For  the 
conservative  case  the  corresponding  equation  is  f(x)  -  u  =  0.  The  reader 
can  easily  verify  that  a  singular  point  (x,0)  is  a  vortex  if  f'(x)  »  0 
and  a  saddle  if  f'(x)  <0.  A  separatrlx  is  a  solution  trajectory  vhlch 
tends  to  a  saddle  point  as  t  approaches  +  »  or  as  t  approaches 

In  the  following  chapters  the  author  will  present  results  for  a 
system  which  describes  the  notion  of  a  soft  spring  (negative  cubic  tern 
in  the  Duffing  equation  ).  The  domain  of  controllability  will  be 
described  and  a  procedure  for  obtaining  the  maximum  number  of  switchings 
will  be  presented.  Also  some  peculiar  behavior  of  the  switching  locus 
will  be  discussed  and  analyzed. 
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CHAPTER  II 

DOMAIN  OP  COWEROLLABILm 


2  Probler.  .fteteacnt 

An  initial  state  of  the  system 

x  +  f(x)  =  u 
or 

*-y 

y  «  -  f  (x)  +  u,  with  -1  <  u  <  +  1  (1) 

shall  be  steered  to  the  origin  in  minimum  tine.  The  synthesis  problem 
is  a  aired  by  determining  the  feedback  control  u(x,y).  The  following 
conditions  ere  placed  on  f(x)» 

(1)  f(x)  Is  in  the  class  (P‘, 

(ii)  -  f(-x)  -  f(x). 

(i*i)  »  fJ(x)  >  0  in  some  region  about  the  origin,  this 

implies  that  there  exists  some  d  >  0  such  that  f(x)  >  0 
for  0  <  x  <  d. 

(iv)  There  exists  cone  x  >  0  such  that  f(x)  <  -  1  for  x  >  x  • 

A  possible  graph  of  f(x)  is  shown  in  figure  1.  Condition  (iv) 

assures  that  each  system  (1+)  with  u  *  +  1  and  (1__)  with  u  «  -  1,  as 

described  in  section  1.3,  has  at  least  one  critical  point,  x  and  x 

8+  - 

a  saddle.  The  origin,  which  la  the  target  for  this  problem, is  a 

stable  e<iullib.'ium  point,  by  condition  (iii),  since  there  u  is  set  to  zero. 


For  this  investigation  the  following  f(x)  will  be  considered: 
f(x)  -  X  -  p  X3  , 

the  charocte ratio  of  a  soft  spring  with  p  >  0  as  parameter. 

System  (1)  can  be  Integrated  in  the  phase  plane  by  eliminating 
the  time: 

(3) 

then 

Tar  u  ■  constant  one  can  exprees  the  trajectories  y(x)  analytically  as 

y(x)  -  +  JyZ( 0)  +  2  u  x  -  2  J*  f(?)  d?  (4) 

0 

With  equations  (3)  and  (4)  and  the  specified  f(x),  x(t)  and  y(t)  can 
ha  determined  analytically.  The  time  needed  to  steer  a  point  (x  ,y(x_)) 

C  C 

to  (x  ,  y(x))  is  by  equation  (3): 

t.J*  s* 

V 

with  f(x)  ■  x  -  p  x3  and  equation  (4) 

t-f  -  ** 


JyZ{ 0)  +2ux  -  x2  +  |  x 


(5) 


which  la  an  elliptic  Integral  of  the  first  kind. 

The  inversion  of  this  elliptic  Integral  [2,  pp.  452-455],  with 
y2(x)  *  h(x),  yields 


1  v. 
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vhare  xf  i?r  "era  of  h{x);  h#(x)  «  .  6(t)  Is  the  Weierstrass 

function.  ‘Its  <)5fiuttior  <  r  d  characteristics  are  presented  in 
[2,  Cu.ptev  /).],  b(o)  -  -  £s;  ~  • 


2.2  1a-]c-.1p  r'r  C2RtrojJ.oTi.lit;  ,  Existence  of  an  Optimal  Control 

As  Mentioned  iu  lection  l.li  in  this  section  the  domain  of 
controllability  for  system  (1)  vith  conditions  (i)....(iv)  an  f(x)  will 
be  described,  end  the  existence  of  an  optimal  control  will  be  proved* 
Theorem  2,?.  Given  the  system 

*  *  y 

y  *  -  i(x)  +  u  (1) 

where  u  ir  rrsouwble  on  -1  <  u  <  +  1  with  f(x)  as  described  in 
ccctj.cn  fi.l  satisfying  conditions  (i)...(iv).  The  domain  G  of 
controllntdlVcy  i»  she  o per.  connected  band  B  between  the  principal  »«“ 
narntricen  through  and  x  .  Furthermore,  each  point  in  B  can  be 

S*  8** 

steered  o  the  origin  by  en  optimal  control  u(t)  6A. 

The  principal seperA trices  arc  the  trajectories  through  x  .  and  x 

0+  G- 

indicated  vs  solid  lines  in  the  phase  plane,  figure  1. 
iTfJi :  Ly  theorem  1,1  there  exists  a  neighborhood  of  the  origin  which 
Uci  in  G,  Thus,  one  au3t  ahev  th-'  _  «....»**«,  xn  a  can  be  steered  into 
this  neighborhood.  Start  at  point  P  (compere  figure  1).  Choosing  u  =  0 
the  trajectory  follows  o  solution  of 
x  =  y 

y  -  -  r(x) 


At  x  <  d,f(x)  >  0,  therefore. 


IT.i'n  y  will  lr  rjf;t:  £.nd  x  vill  dec reuse, 
to  centime  in  t>  e  azrs-  direction  act  u  -  -f  3..  If  r(;;,  >  1  for  core 
value  of  x  <  d  (nee  cinched  line  for  f(x)  JLr  figure  1)  v  vi31  cccreaco . 

Bo*  the  trajectory  cannot  leave  the  region  P  fcounde  by  ti  e  principal 
nepnratri’:  through  xr+,  that  in  for  r>  ncxin^i  conntont  u  =  +  *..  If  the 
trojcctorj  tent’s  to  intersect  the  positive  x-lxIg  set  u  »  0  anti  y  will 
decrease,  vhi3jC  x  in  3till  decreasing.  ’Phus,  one  con  steer  the  t rejects 
into  the  thii’d  ourlrart  ^tith  a  properly'  chosen  u.  In  the  tail’d  trod  - 
cry  value  0  <  i  <  1  will  steer  tlie  trajectory  tovardr  the  negative  x- 
nxLs.  A3  it  enters  the  second  quadrant  x  and  y  vill  increase.  To  steer 
it  tovord  the  origin  u  =  -  1  ic  needed.  Ac  it  encircles  the  origin  ‘,2th 
decreasing  distance  fron  the  origin  one  reeds  only  to  pick  the  proper 
u  to  conti r re  In.  the  vented  direction. 

If  P  lice  ir  one  of  the  other  quadxxmtc  th'’  proof  is  eirii3.tr. 

ihe  .ealor  should  cl30  observe  that  a  trajectory  *>hich  enters  the 
region  bounded  by  tie  four  sopars trices  closest  to  the  origin 
(ace  curvilinear  qpadrilate-nJ.  inside  the  solid  one.  dashed  c«paratrlc  ■  s 
in  figure  3  )  vill  not  leave  it  again. 

Thu:-,  it  her  been  eh  ova  that  each  initiu}.  point  P  P  3  can  be  rteoved 
to  the  origin  some  control  u(t)  Since  x  end  y  remain  bounded 

during  this  nrooenc  by  theorem  1.2  chore  exists  un  optimal  control. 

This  ccxplcton  the  nr oof. 

Tbo  r.rco^  in  boslcly  the  some  for  the  enac  "./hen  positive  dating 
1  i,  present. 

By  the  rj«xi.mxr  principle  and  the  Theorem  of  Interlacing  Zeros  the 
optimal  controller  is  a  relay  controller:  u  assumes  the  values  -1  and  +3 
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2 . 3  Pha^  flan^  Aaulyris 

Since  the  'behavior  of  the  trajectory  through  ■’.he  origin  in  import;  nt 
for  the  switching  locus  a  detailed  discussion  of  the  -clnticn  ‘between 
*(jc)  and  the  trajectory  through  the  origin  in  pteccrtcd. 

Ptor  this  trajectory  in  equation  (H)  y(o)  -  0  and 

y(x)  *  +  u  x  -  2  J  f(?)  dr 

Vith  the  function  as  described  in  section  2.1  one  rupc  distinguish  thr- 

different  cases. 

(a)  Each  system  (l+)  vith  u  *  +  1  and  (l_)  vlth  u  =  -  1  hts3  one 
oingular  point  as  the  only  reel  solution  of  each  equation 

f(x)  -  l-xa+ 

and  f(x)  — l-xg_j  (4) 

see  also  the  graph  in  figure  2.  The  trajectory  through  the 

origin  (solid  line  for  u  »=  +  1,  dashed  line  for  u  =  -  l)  has 

negative  slope  throughout  the  second  and  fourth  quadrants. 

4 

This  case  applies  to  the  above  example  for  0  >  — r;  - 

—  I 

(b)  Each  system  haG  three  singular  rcirts,  i.e.,  each  equation  c? 
(U)  has  three  real  roots.  The  graph  of  f’x)  and  the  phase 
plane  are  r.hovn  in  figure  3.  For  u  =  +  I  the  seporatriccc  arc 
trajectory  arc  sketched  ir  a  oil’  liner,  for  u  *  -  1  th^y  arc: 
indicated  in  drshed  linen.  ‘Jvo  a-slltst  root,  xQ,  ic  a 
vortex,  vhcruar.  and  xg  rupi'OGent  saddles.  H\e  separatrix 
through  Xj.  encirclea  xft.  It^  rcintereection  vith  the  y-ax’c 
lies  in  the  same  half  plane  as  Aercforc,  the  tmject^iy 
thrjugh  the  origin  dots  not  reintcrcect  the  x-axis.  However 
it  changes  its  slope.  The  slope  is  negative  or.  (-*>,  Xj..), 


FIGURE  2.  ILLUSTRATION  OF  CASE  <o). 
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(x^  xq+),  and  ®nd  positive  on  (x^._,xo_  )  and 

at  x^  and  x^+  the  slope  is  zero. 

2  4 

This  holds  for  the  example  when  —  <  3  <  7^.  The  lim'  ting  case 
between  cases  (0)  and  (b)  is  the  case  when  f(x)  Juec  touches  the  lines 
u  *  +  1  and  u  *  -  1.  Each  equation  (4)  has  three  real  roots,  a  simple 


root  x  and  a  doable  root  x  * 


When  x  has  the  value  of  this  double 


s - o  V 

root  the  trajectory  through  the  origin  baa  2ero  slope,  except  at  these 

two  points  it  has  negative  slope  throughout  the  second  and  fourth  quadrants . 

4 


For  the  example  0  - 


27* 


(c)  Each  system  has  three  singular  points  aa  in  (b).  But  the 
aeparatrix  through  reintersects  the  x-axls  in  the  other 
half  plane,  see  figure  4.  So  the  trajectory  through  the  origin 
lies  Inside  the  region  surrounded  by  the  sepertrix  through 
X.  and  therefore  reintersects  the  x-axls*  This  holds  for 

?<!r 

The  limiting  case  between  cases  (b)  and  (c)  is  the  case  when  the 
ceperatrix  through  reintersects  the  x-axis  in  the  origin,  that  is,  the 
trajectory  through  the  origin  coincides  with  this  aeparatrix  and 

2 

reintersects  the  j^axis  in  a  ingular  point;  this  occurs  for  p  » 


1 
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CHAPTER  III 

SWITCHING  LOCUS  FOR  SMALL  VALUES  OF  P 


In  this  chapter  the  author  vlll  dlecuao  the  cnee  (c),  l.e.,  for  the 
_ rv  only  values  (S  <  vlll  he  considered.  Oen  for  0  <  p  <  §f 

one  can  find 

1  <  xQ  <  1»1 

Once  xQ  ia  date  mined  one  obtains 


for  mc.il  p  one  can  approximate 


‘o  •  x»  v-Jj  ■  i  ■ 

3.1  LqcuB-M*^  °f  Switchings 

por  the  given  system  (1)  the  theorem  of  interlacing  ceros  which 
y„.  ateted  In  aectlon  1.4  as  basic  for  the  construction  of  the  switching 
loeu.  is  equivalent  to  Stum’s  caparison  theorem  for  dlfferentlel 
equations  [4J.  for  one  «s  write  the  second-order  differential  equations 

for  y  and  t2: 

j  +  f'(x)  y  *  o 

^  +  \  *  0# 


and 


HI i  —  er*. 


As  stated  in  section  lA  the  segments  of  the  trajectory  through 
the  origin  in  the  second  end  fourth  quadrants  belong  to  the  switching 
locus ,  which  ere  and  W^.  One  follows  the  procedure  described  in 
lA  to  obtain  consecutive  arcs  of  Burmeister  presents  a  different 

procedure  for  the  construction  of  the  switching  locus  in  [5l»  for  which 
he  applies  the  variation  equations  of  the  system.  He  claims  that  his 
method  is  quicker  and  more  accurate  since  fewer  integrations  are  performed. 

The  author  compared  Burmeister ’  s  method  with  the  one  described  in  lA  for 

•  •  ^ 

x  +  x  +  2x  =  u.  The  former  yielded  only  two  arcs  of  the  switcning  locus 

of  desirable  accuracy  and  took  relatively  more  computer  time  than  the 

latter  which  yielded  four  arcs  of  desirable  accuracy. 

If  a  switching  lies  on  the  x-axic  (a  zero  of  TL,  is  a  zero  of  y) 

then  the  next  switching  must  also  lie  on  the  x-axis.  Thus,  consider  the 

arc  in  the  fourth  quadrant,  it  re  intersects  the  positive  x-axis  in 

x^+,  figure  5.  Starting  at  this  point  backwards  in  time  with  u  »  -  1 

the  trajectory  will  reintersect  the  negative  x-oxis  in  x^_  which  1g  also 

the  endpoint  of  .  The  trajectory  leaving  x^m  with  u  »  +  1  will 

■a 

reintersect  the  positive  x-axis  in  x.^+,  the  endpoint  of  W^,  provided  the 
trajectory  lies  still  inside  the  region  bounded  by  the  seperatrix  through 
(see  dashed  line  in  figure  5).  If  it  lies  outside  this  region  (in 
figure  5  the  trajectories  leaving  x^_  and  x^+)  it  will  stay  between  the 
principal  seperat rices  and  go  to  infinity  in  G. 

This  shows  a  way  of  determining  the  maximum  possible  number  of 
switchings,  namely,  follow  the  trajectories  whose  switching  points  lie 
on  the  x-axis  until  they  leave  towards  infinity  and  count  the  intersections 
Xj,  on  the  positive  or  negative  x-axis.  This  number  of  intersections 
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t 


(including  the  origin)  is  then  the  maximum  number  of  switchings.  For 
the  above  example  f(x)  c  x  •  p  x  with  p.  <  —r  the  following  algorithm 

2  I 

has  been  developed: 

The  trajectory  through  the  origin  for  u  *  +  1,  W^,  is 
y(x)  *  -  J |  x*  -  x2  +  2x 

Its  reintersection  with  the  x-axis,  x^+  =  x^  is  determined  by  solving 
^  -  x2  +  2x  e  0, 

that  is  finding  the  smallest  positive  root  of 
^x^-x  +  2  *0  . 

Ct 

Mote  that  by  the  symmetry  of  f(x)  Xj_  *=  -  x^.  lhen  the  next  step  is 
to  determine  the  intersection  Xg  for  the  trajectory  leaving  -  x^. 

This  trajectory  is  y(x)  «  -  /  ^  x^  -  x2  +  2x  +  4x^,  and  Xg  is  determined 

by  solving 

|  x*  -  x2  +  2x  +  4xj  =  0 

for  the  smallest  positive  root.  Then  the  n-th  reintersection  is 

found  with  the  knowledge  of  the  preceding  intersections  by  solving 

.  n-1 

^  x4  -  x2  +  2x  +  4  £  x.  =  0  . 

*  i=l 


Since  x  ,  is  known  and  -x  ,  is  also  a  root  of  the  above  equation  ■'■his 
a-i  m 

fourth  degree  equation  can  be  reduced  to  a  cubic  equation 


*3  *  Vl  *2  -  (t  -  x„-l)  x  +  (5  +  f  xn-l  -  x„-l3)  ■  °» 

which  can  be  solved  exactly;  x  is  the  smallest  positive  root.  However, 

n 

for  computer  solution  the  above  fourth  degree  polynomial  can  be  solved 
more  quickly  by  Newton’s  method  (or  any  other  algorithm)  with  a  starting 


velue  =  >:n_i  +  ^  than  the  cubic  equation  vhich  -*nvolves  variable 
transformations.  And  one  also  must  determine  the  smallest  root  of  the 
three  roots. 


The  above  procedure  is  terminated  when  x approaches  when  the 
trajectory  leaving  -  x  ^  reintersects  the  positive  x-axis  between  the 
reintersection  of  the  seperatrix  through  x^.  (u  -  -  l)  and  x^+.  For  the 
trajectory  leaving  the  corresponding  point  »x  the  polynomial  has  no 
positive  root  x  <  x^ 

These  trajectories  with  the  switching  points  on  the  x-axis  divide 
the  domain  G  of  controllability  into  subregions  for  each  of  which  a 
certain  number  of  switchings  holds,  e.g.,  in  figure  5  for  each  point  in 
a  region  marked  with  n  *=  4  four  switchings  are  necessary  to  steer  it  to 
the  origin  in  minimum  time.  The  outer  boundary  of  each  region  belongs 
to  it,  except  the  boundary  of  G,the  seperatrices  through  xg+  and  xg_, 
which  does  not  belong  to  G.  The  switching  at  the  origin  to  u  =  0  is  not 
included.  The  example  shown  in  figure  5  is  plotted  for  0  =  0.01. 

For  the  linear  case,  0=0,  the  harmonic  oscillator,  the  difference 

between  two  succeeding  intersections  of  the  switching  locus  is 

x,  -  x,  ,  =  2  and  x  =  2n.  For  0  >  0  it  is  x.  -  x.  ,  =  2  +  6.  with  5. 
i  1-1  n  I  l-x  l  i 

increasing  for  increasing  i.  In  Markus -Lee’s  example  which  corresponds 


to  0  <  0  the  difference  x^^  -  x^_^  is  decreasing  for  increasing  i.  As  a 

2 

rough  estimate  for  the  maximum  number  of  switchings  for  0  <  0  <  ~  one 


may  assume  the  linear  case  and  check  how  many  intersections  with 
xi  “  xi-l  =  2  are 


possible  on  (0,jO  or  n  <  —  (for  very  small  0 


n<^Vr 


The  switching  locus  of  an  example  with  0  =  0.03  is  shown  in  figure  6. 


Here  the  maximum  number  of  switchings  Is  two 


' 


In  figure  7  the  Isochrones  are  shown  for  the  preceding  example 


(0  ■  0.03).  Their  equation  is 

*(*.*>  -  f  -  constant  ,  for  the  tine-optimal  u. 

Sadi  point  on  each  a  dosed  curve  is  steered  to  the  origin  tine-optlnally 
in  the  sane  length  of  time. 

These  Isochrones  represent  also  boundaries  for  the  sets  of 
attainability  with  the  origin  as  initial  point,  e.g. ,  the  isochrone 
with  T(x,y)  -  Tk  surrounds  a  region  each  point  of  which  can  be  reached 
froa  the  origin  with  the  tine  optimal  control  in  the  tine  t  <  T^. 


from  this  illustration  and  other  examples  one  nay  conjecture  that  these 
regions  are  convex  if  they  lie  within  x^._  <  x  <  x^+. 


3*2  Influenoe  of  Dating 

The  presence  of  damping  affects  the  behavior  of  the  systen 
significantly.  Since  the  nonlinear  systen  with  damping  cannot  lie 
analysed  the  linear  systen  will  be  considered  first.  The  differential 
equation  satisfied  by  a  tine -optimal  control  is  given  by*i  +  +  x  =  sgn 

If  x  ■  y  then,  Y  +  of  *  y  -  0,  and  the  equation  satisfied  by  tj2  is 

*2  '  °  '  0 

The  systen  oscillates  only  for  a  <  2  (a  *  0  represents  the  harmonic 
oscillator)  so  for  a  >  2  the  nsrlnon  maker  of  switchings  is  one, 
because 


backwards  in  tine,  with  the  only  sero  at  t  ■  0.  For  a  ■  2 

a  t 

T)g(t)  -  (At  ♦  B)  e  , 


where 


IMs  number  may  be  changed  by  1,  depending  on  whether  the  point  with  xn 
as  abscissa  lies  above  or  below  the  switching  locus. 

With  the  nonlinear  term  the  equations  for  y  and  ^  are 

y  +  cr  y  +  f'(x)  y  =  0 

\  '  7  ^2  +  f'M  V2  *  0 

with  f'(x)  =  1  -  3Q  x2 

2 

The  solution  oscillates  for  f  '(x)  -  p-  >0  or  for 


A  rough  estimate  on  the  maximum  number  of  switchings  can  be  obtained 
by  substituting 


in  the  above  inequality  for  n. 

In  figures  8,  9 >  and  10  an  example  for  3  -  0.005  is  shown  to 
demonstrate  how  the  damping  changes  the  shape  of  the  switching  locus 
and  the  maximum  number  of  switchings. 


3  t  ,  * 

With  Z  =  e  y  one  obtains  for  the  above  equation  Z  +  (f  (x(t))-  qp)Z“0 

for  which  the  comparison  theorem  [4]  can  be  applied. 


FIGURE  9.  SWITCHING  LOCUS  FOR  £*0.005,  WITH  DAMPING  a  *0.2. 
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CHAPTER  IV 

SWITCHING  LOCUS  FOR  INTERMEDIATE 
AND  LARGE  VALUES  OF  0 

In  this  chapter  the  cues  (a)  and  (b)  as  described  In  section  2.3 

will  be  discussed.  Originally  the  aasusption  was  aade  that  If  the 

trajectory  through  the  origin  does  not  relntersect  the  x-axls  then  It 

oelnoldes  vith  the  switching  locus  over  Its  entire  length,  and  the 

miier  of  switchings  is  one.  Bat  computer  results  showed  that 
2 

there  are  cases  for  0  >  ^  for  which  there  exists  an  arc  of  the  switching 
locus  sticking  back  firoa  infinity.  Furthtrnore  the  switching  locus  Is 
41s cent lnoous,  the  arc  through  the  origin  Is  not  connected  with  the 
eae  sticking  back  from  infinity.  Figures  11,  12  and  13  show  exanpi.es . 

ISo  different  integration  procedures  were  applied  to  asstaae  the  result. 
Oae  integration  was  with  respect  to  tine  t,  the  other  procedure  was  with 
respect  to  x  as  discussed  later.  Froa  these  results  the  conclusion  was 
draws  that  there  exist  arcs  of  the  switching  locus  sticking  back  froa 

2  ti 

infinity  only  for  sane  subinterval  of  —  <  P  •  In 

the  fallowing  sections  the  author  will  consider  certain  situations  for 
2 

0  >  —  and  show  whether  there  exists  a  second  switching  or  not. 
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b.l  Derivations 

The  systea  equations  are 


i  -  y 

f  -  -  f  00  +  u 


(1) 
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FIGURE  12.  SWITCHING  LOCUS  FOR  /9*0.08. 
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vlth  u  *  +  1  below  an &  u  ■  -  1  ibovt  the  twitching  loco*.  The  adjoint 

lyitia  la 

\  -  f'to  1L 

(2) 

and.  the  Haadltonian 

-  \  y  +  \  fn  “  f00] a  h  »  (6) 

share  h  it  a  pot  it  ire  constant  along  an  optlaal  tratf  ectory.  The  equation 
of  an  optlaal  trajectory  vat  derired  in  lection  2.1  at 


y(x)  -  +  Jy2(0)  +  2  u  x-  2  f  f(5)  d?  . 


00 


(2') 


The  tana  procedure  la  applied  to  obtain  TV^  and  7^  at  functions  of  x  and 
J 

S’  V"y  ^ 

d1>2  \ 

S’  \  ’  *  T  • 

One  usee  the  Haslltonlan,  equation  (6),  to  ellnlnate  7^: 

With  u  -  f(x)  »  yy'  froa  tycten  (l)  one  obtains  a  first-order 
differential  equation  for  ^  In  x. 

*  ^  -  *'  \  ’  -  I 


which  can  be  Integrated 


(7) 


•  3^* 


Vxo> 

S  W' 


1  -»£*-  ,  y  /  °* 

Jo  y3(?) 


Thua,  one  can  write  the  feedback  control  in  analytic  form 

*  ngn  (y(x)  [^4  -  »  f-f-}  >>*°- 


n,e  eoostante  xo  end  h  are  specified  in  the  next  section.  With 
f(x)  .  x  -  0  x3  the  above  Integra  is  an  hypewlliptlc  Integra  of  the 
first  an  [3]  which  cennot  he  solved  in  closed  form. 


4.2  Switching  Condition 

*o  show  the  existence  or  non-existence  of  *  second  switching 
ODe  Mt  investigate  station  (7).  ">1.  is  done  the  following  VW. 

The  lnitia  conditions  are  chosen  on  the  switching  locos  in  the  second 

^tant  fw  which  o  -  -  1,  •«  «««•  l1*-  *»  u 

patorm 4  Oong  the  trajectory  hOnw  the  switching  locos.  B>»s,  **» 

now  on  only  the  case  o  -  ♦  1  I*  considered.  Sr  the  sjr—etry  of  f(x) 

the  eitoation  is  exotl,  the  ...  Shore  the  etching  loco,  with  opposite 
ago.  let  x  «  -  a,  then  the  lnitia  ordinate  in 


jK-a)  e  ,f  2a  - 


T 

a 


Since  y(-n))  is  on  the  switching  locos  hg(-a)  -  °”  ““  choo,e 

y_e)  erhitrarily  hocouae  syeten  (2)  is  hooogeneous.  mtt  he 

~4d.ee>  Orm  -1  to  +1  at  this  initial  point.  Thus 
positive  far  u  switches  m*  -a  °  ^ 


lat  \(-»)  *  1  th*n  **  e<IUAtlon  W 
h  *  y(-e) 

There  exists  a  second  switching  when  the  integral  in  equation  (7) 
become  xero  for  sem  x  ^  -a. 


Homver,  Uo»«n  see  fro.  figure  rW  i»  ** 

„  that  the  int.gr.tion  IntervU  »t  he  divided  into  «W»terval.  for 
rtleh  y<x)  1*  single-valued.  The  first  interval  i.  t-M.1  fOT  ’’hl<* 
y(r)  H  in  the  .econd  quadrant,  .here  re  is  the  ahociie.  for  .hich 
y(r  )  *  0.  The  second  interval  is  (xc,x3  f°r  vnich  y(x)  is  in  h 
Imrer  >'*1  r  plane.  The  required  division  hy  y(xc)  causes  dlfficu 
ut„  on  in  the  evaluation  of  the  Integra  it  *1U  he  shovn  ho.  to 

avoid  this  problem 

from  (6)  one  can  determine  1^(xc) 

h  -  \{*c)  y(*c) +  Vxc>  ^  "  f(xc>l  * 

Since  y(xe)  H  0 


Vxc>“I^T  » 


T^(xe)  >  0  since  1  -  f(*c)  >  0  tor  x|+  <  x£  <  0. 

&  _  ...  aw.  __vt  witdxinx  -  if  there  exists  one* 

It  will  be  shown  Inter  that  the  next  ewxtcaine 

occurs  in  the  fourth  quadrant. 

Therefore*  first  nil!  he  dotandnad,  sad  than  the  «matlon. 


V*>  „  .  b  { 

-y^xT  H5T 


x  >  0 


nill  be  exaadLned.  At  a  second  switching  point  ^(n)  *  °> 

^.hf  -|5--o 

J0  P(?) 

.01^  for  x  ^11  glv.  a  point  (*,y)  of  th.  etching  loeus.  In  B.neral 
thla  egaatlon  cannot  ha  solved.  But  if  there  is  a  sidtehlug  for  som. 
o  <  X  <  -  the  expression  on  ths  left  hand  side  of  this  equation  .ill 
Hup  it.  sign.  Since  Tt,(x)  >  0  and  y(x)  <  0  in  th.  third  quadrant, 


y°) 

“Hoy 


<  0.  Thus,  there  exists  a  switching  on  (0,»)  if  the  expression 


TL(°) 


o  y^(x ) 


is  positive,  or  if 

t  *  dx 


*’  o  yJ(x) 


y°) 

<1TW1  • 


Considering  absolute  values  of  jy  the  switching  condition  is  given 

A 


r  dx  .  2v  ; 

Jol7wT  >*mr. 


V°) 


f  dx  .  2'  ' 

”  J.pwT  ^ 

there  Is  no  switching. 

4.3  Estimation 


In  this  section  estinates  are  derived  to  examine  the  inequalities 

y°) 

(9)  and  (10).  is  determined  by 


V°) 


<n  (x)  0 

u.  r  -£!  -  h  f  -js- 1 . 

*<  -  'xy3(?)  J 


where  the  integration  is  performed  in  the  third  quadrant.  The  above 
integral  is  subdivided  into  two  integrals  as  follows 

f°  **  ■=  f*&  dg  +  f 

x  y^(V)  x  y'^CO  -ay  (?) 

Purthermore  with  y  y'  =  1  ■  f(x)  the  first  integral  of  the  right  hand 
side  is  transformed  as  follows 


r-s--  f 


(-a)  dy(?) 


x  y3(?)  Jy(x)  y^(5)D--«*(?)3 
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Bnu,  the  equation  to  be  investigated  is  cow 

V*>  m  rVl*  h  f(~*^ _ i.v,  r°4s- 

rT5I  '  ^  LH*T  Jy{l)  7^7-1 J  J-»  y3^)  ■ 

0 

for  the  Mtimtion  of  1  -  f(x)  on#  need*  to  distinguish  three 
different  esses,  ss  is  showa  in  figure  15*  At  x  »  -  **** 

its  audJHt  Tslne. 

Csse  It  **f<xe<““<“a<0 

^  c  /53 


1  -  f  (x)  >  1  -  f(-n)  on  (-  “»"•) 

1  -  f(x)  >  1  -  f<ic)  « (v  •  ^ ) 

This  refKlres  thet  y(-  — )  exists: 


w*- 

,  the  condition  on  s  is 


iBe  /3e 


*(-2-  ♦  -*-)<. <-i- 

*  'ifla  7  /35 


**(-£-  +  ^.')  +  JC(-5 - §-WthO<y<l. 

*  v»  ifle  7  *  73*  Oflp  7 


2* 

•+  c 


then  1  -  f(*)  >  1  -  *(*e)  ••  (xe> 

Therefore,  s  nut  be  sndj  thnt  )ze  j  <  jz#  |  exists.  Bms 

JL  <  s  t  |*g|  +  33  . 

J3 
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'wl»* 


Case  3: 


/36 


<  x  <  -a  <  0 
c 


then  1  -  f(x)  <  1  -  f(xc)  on  (*c,  -a). 


Tne  liinita  on  v.  are  then 


0  <  a  <  t  (—  -  ~— 

*  Vis  183  ' 


For  case  1  the  integration  ic  performed  over  the  intervals 

—) 

/33 


i  y(x),  y(-  |  and  I  y( - -)  ,  y(-a)  I  .  Then  the  estimate  is 

L  /B  J  1 


V0^  ,,  fVx^  h  r  i  in 

7T5T  <  L‘x  Vy^T  *  i^T  OT7  "  “ T“  1 

c  y(-~) 

/3s 


X*+X_ 


.  0 


.  -*■  f— - L_  .  h  I  . 

1“f^  ^  ^y(-  —  )  y(-a)  J  -a  y3(x) 

/To 

After  performing  the  limit  and  taking  absolute  values  Ter  y 

¥°>  .  r 0  te  ,  W 

WI  ha  y^(x)  h{l-f(x„)1 


a  -  r(*cn  !y  )l 

c  /3s 


i  -"f'TO 


[ 


|y(-  —  )| 

/36 


-k) 


(13) 


This  holds  for  ease  1.  The  estimates  for  cases  2  and  3  are  given  later. 

In  the  above  inequality  (13)  the  integral  and  T\L(XC)  n3UJS^  be 
estimated.  The  equation  for  T^(x)  can  easily  be  found  from  equation 
(7)  and  equation  (6). 


i 


»<S?  -• 


Startii*  «(t  xo  -  -  a  on  the  switching  locus  the  integration  is  performed 
on  the  trajectory  in  the  second  quadrant.  Fallowing  the  above  estinatlon 

yrooodure  one  obtains  for  case  1 


far  case  2  T^(xc)  >  1 

I 

for  eaae  3  \<*0)  < 

m  the  fallowing  estiaatea  are  pursued  which  Show  that  there  is 
no  eessnd  switching,  by  inequality  (10),  for  cases  1  and  2. 

For  the  integral  in  (13) 


the  cubic  function  -  2*  4  is  replaced  by  its  tangent  at  x  »  0 

)4 

which  is  4-2x.  Then  the  integral  can  be  estlnated 


o 


Finally,  the  desired  estimate  for  case  1  is 


The  estimate  ‘or  case  3  is  given  on  page  ;9. 

How  one  needs  an  estimate  for  the  1-  ft  side  or  inequality  (10). 


V'x) ' 


o  ’■*  (°>xJ  -  ?.y.  -ti)  tin 

l  c  J 


•Hie  cubic  polynomial  is  approxinm  od  by  i‘  r.  e anger  .s  ir  ii-lust rated  ir 
figure  IS.  Then 


-1*7- 


After  perforaing  these  integrations  one  obtains  the  folxow  ng 
inequality  as  equivalent  to  the  condition  that  no  second  switching 

exists  in  case  1. 


[  r  -Wt  - 1  1)  •= 0 


Th.  ecconi  braces  contain  the  quantity  obtained  In  Inequality  (it).  It 
b.  replaced  by  the  one  from  Inequality  (15)  if  =»se  2  i8  examined. 
All  square  roots  In  the  first  braces  exist  for  B  >  .  Th*  function 

-  f(a,0)  is  plotted  versus  0  in  figure  IT  vith  y  nc  paraneter 

lies  in  the  shaded  area  for  0.2  <  y  <  1.0  and  B  >  ^  •  This 
shorn  that  Inequality  (16)  bolds,  l.e.,  that  there  is  no  second 
switching  for  9  >  end  y  (jig  *  <  “  /5@  ' 

For  case  2  a  positive  tern  in  ?(a,B)  drops  out,  so  that  the 


gpproodmtion  becooas  even  better. 


-*9- 


g 

To  show  the  existence  of  a  second  switching  for  some  0  >  ^ 
and  none  a  inequality  (9)  nuat  be  verified.  Since  there  is  apparently 
ealy  a  snail  range  of  a  and  0  for  which  a  second  switching  nay  exist 
(sae  figures  12  and  13),  and  since  the  estimates  are  not  very  sharp  one 
example  will  be  tested,  which  will  be  a  -  0.2,  0  -  0.06,  see  figure  12. 
One  can  see  that  case  3  applies:  -  <  x  <  -a,  for  which 


V°)  _  2  f 0  di 

ktttsti  »  1. ^  ■ 

O 

There  the  cubic  function  ^-2x  +  0x  is  replaced  by  its  segments  over 
8  intervals  on  -a  <  x  <  10.  After  the  evaluation  of  all  the  estimates 

one  obtains 


dx 

|y3W  i 


>2.9 


V°>  <2ll 

* 

satisfying  inequality  (9).  Thus  there  exists  a  second  switching  point. 
Because  dr}2/dx  exists  and  xs  nonzero  if  r?2  =  0  there  nust  be  an  interval 
about  a  *  0.2  for  each  point  of  which  a  second  switching  point  exists, 
lhe  sane  can  be  said  for  0. 

<4  Bstlaation  for  Snail  a 

To  prove  the  results  in  figures  12  and  13  one  wist  show  that 
thare  is  no  second  switching  for  snail  a.  A  different  approach  is 
used  to  verify  inequality  (10).  x  is  expressed  as  a  function  of  y, 
and  the  integration  is  psrforasd  with  respect  to  y.  The  requirement 
of  a  being  —yi?  implies  that  the  range  of  x  and  y  is  smll  in  the 
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left  helf  phs.se  plane,  for  the  equation  describing  the  trajectory, 

2  .  _  2  ft  4 

y  *  4a  +  2x  -  x  +  ej  x  > 

one  can  apply  an  approximation  procedure  to  obtain  x  as  a  function  o 

y 

i — i — nr 


/  2  2  B  4 

-  7 r  ■ y  +  *  Vi 


2 


where  r  »  1  4  4a.  This  sequence  converges  for 

8  4  2  2 

£  x  <  t  -  y 


This  approximation  results  in 

2  U  6 

x  *  x  +  A  y  +  By  +  C  y  +•  •  • 


4  - 


32r' 


r  r' 

a*  /a 

£j 

> 

r~  r 

r 


B  “ p  -  f  ^ 1 (i^-*  •  ••• 


+  - 


r 

4  .2 


and  xc  «  (1-r)  -  g  +  §5 


.8 


♦  -  #  < 


2  3 

-2a(l  -  a  +  2a  +  20  a  . 


The  higher  terms  have  coefficients  with  decreasing  aagnitude.  The 
above  expressions  are  evaluated  by  developing  series  in  a  since 

,  whose  power  series  converges  for  4  &  <1.  The  following 
estimation  is  only  valid  for  sufficiently  small  a.  The  integration 
in  the  second  quadrant  in  then 
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U*) 

j(x) 


bj  JS- 

-»  y3(?) 


r  xc  *  *  xc'x  ,  ,y(x)  Vs  „.i 

■ ' [  Pw  ’  7w  ‘ 3  J*.,  7"  ** 


After  substituting  the  polynomials  in  y  and  performing  the  Integration 

obs  obtains 

4ik  ?  k  ? 

"  yT^7  "  ^  +  UBh'  "  2CH  -  ^Bhy(x)  -  2Chy*(x). 

Lotting  x  approach 

11m  [Hg(x)  -  2Ah]  *  -2A  +  Ugh2  -  2Ch\  (17) 

X-r»xt 


Thu, 


3  4 

T^g(xc)  =  2Ah  =  h  «  0  ^  i  which  compares 


with  Tk,(x)  =  "j  in  the  first  three  terns  of  its  series  in  a. 


In  the  third  quadrant  the  integration  is 


V°>  Vx>  .  f  d? 
^r  =  H^r'hJx^ 


H7T 


-  h 


-  3  f (0)  v!  ^ . 

"y3(x)  y3(0)  Jy(x)  y 


With  the  abort  polynomials  and  equation  (17)  one  obtains  the  following 

estimate 

sTyrsrr  *  r  f1  *  tu]  • VB  Ai  [x  *  • 2C{k*? ri  4  “N/sl4-" 


(■*.) 
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or  as  a  function  of  a 


yo>  1  r 

KH5TT  -  75  L1 


5.29a  +  10 


?. 

a  + 


For  the  integral  in  the  fourth  quadrant  ore  can  use  the  estimate  for 
case  1.  Then  the  condition  that  there  is  no  second  switching  for 

small  a  is 


[ 


■f*  -  2^fs  •  "§7b  -h  *  7-3(^B'  • 


-3 


°-8’  Jm*1  -1 

*  * 

<  (l  -  2.46a  -1.3a2  +  .  (18) 

Tils  inequality  is  illustrated  in  figure  18.  The  left  hand  side 

of  inequality  (18)  is  plotted  versus  a  with  0  as  parameter.  The  right 

hand  side  marks  the  boundary  of  the  region  for  each  pair  (a,0)  of  which 

the  inequality  holds.  The  inequality  was  also  tested  for  3  =  0.09  and 

was  satisfied  for  o  <  0.1.  In  the  preceding  section  it  was  shown 

that  there  is  no  second  switching  for  a  >  ^  )•  For  9  >  2 

(18)  holds  for  all  a  less  than  the  above  value.  This  inplies  that 

there  is  no  second  switching  for  0  >  2.  There  does  not  exist  a 

2 

second  switching  for  any  0  >  ^  and  sufficiently  small  a. 


4a  -  (1  - 


72 76 


?£l- 
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4 . 5  Non-Existence  cf  Switchings  In  the  Third  ffltadrant 

The  trajectory  leaving  the  switching  locus  through  the  origin  in 
the  second  quadrant  backward  in  time  is  given  by 

y+(x)  =  +  J* a  +  2x  -  x  +  |  x 

Assume  TV,  =  0  in  the  third  quadrant  at  x^.  Then  u  switches  from 
+1  to  -1  and  the  trajectory  is  now  described  by 

y.(x)  -  +  ^(e  4  xx)  -  2x  -  x"*  +lfx*  . 

If  x^  <  -  a,y  (x)  reintersects  the  negative  x-axis  entering  the  second 
quadrant  again,  since  it  cannot  intersect  the  y-axis:  y_(0)  *  J*  +  x1) 
is  imaginary.  Thus,  this  is  not  a  time-optimal  process.  Hie  trajectory 
is  now  closer  to  the  origin  than  when  it  started  going  backwards.  If 
-a  <  the  trajectory  y_(x)  intersects  the  y-axis  at 

7_(0)  =  -  +  xj  >  -  Jit  , 

where  -  =  y+(0).  As  the  trajectory  y_(x)  enters  the  fourth 

quadrant  continuously  increasing  it  must  intersect  the  switching  locus 
throug  i  the  origin  in  this  quadrant  at  x2  , 

x2  =  a  +  x1  . 

with  -a  <  x1  <  0:  0  <  xg  <  a. 

A  trajectory  corcing  from  above  this  switching  locus  in  forward  time 
reaches  the  origin  more  quickly  along  this  switching  locus  than 
along  the  prescribed  way.  It  takes  evei'  longer  to  reach  the  origin 
from  the  starting  point  (-a,  y(-a))  than  from  (x2,-y(x2)).  Thus, 
there  cannot  be  a  switching  in  the  third  quadrant  for  a  time-optimal 


control. 


CHAPTER  V 


SUMMARY  AND  RECOMMENDATIONS 


In  this  investigation  results  for  the  synthesis  problem  were 


obtained.  It  was  shown  that  the  coefficient  3  >  0  of  the  negative 


nonlinear  tern  of  the  Duffing  equation  characterizes  the  behavior 


of  the  system 


the  maximum  number  of  switchings  kis<»>k>2 


A  procedure  was  developed  to  divide  the  domain  of  controllability  G 


into  subregions  for  each  point  of  which  a  certain  number  o 


switchings  is  necessary  to  steer  it  to  the  origin  time-optimally. 

With  this  procedure  one  also  obtains  the  maximum  number  of  switchings 


Computer  results  opposed  the  assumption  that  for  3  > 


trajectory  through  the  origin  i.i  the  second  and  fourth  quadrants. 


coincides  with  the  switching  locus  over  its  entire  length,  implying 


that  a  *  1.  Ntaiely,  arcs  of  the  switching  locus  sticking  back  from 


Thus  a  switching 


infinity  were  discovered  for  values  p  > 


condition  involving  functions  of  x  and  y  was  derived.  Since  it  contains 


a  hyperelliptlc  Integral  it  was  necessary  to  derive  estimates  to 


examine  this  condition.  Using  these  estimates  it  was  shown  that  there 


exists  a  second  switching  for  an  Interval  about  3  =  0.09  and  about 


a  =  0.2.  a  is  the  abscissa  of  the  point  where  the  trajectory  intersects 


the  switching  locus  through  the  origin.  It  was  also  shown  that 
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there  does  not  exist  a  second  switching  for  sufficiently  snail  a 
2  4 

and  any  9  >-%z  ,  and  if  fl  >  —  there  does  not  exist  a  second  switching 

M  W 


5.2  Recommendations 

With  these  estimates  one  could  rot  -'over  the  whole  range  of  a 
to  determine  whether  there  exists  a  second  switching  or  not  for 

4 

•gj  <  0  <  2.  Thus  more  investigation  on  this  hyperelliptic  integral, 
equation  (7),  is  necessary  to  obtain  more  information  about  the  arc 
of  the  switching  locus  sticking  back  from  infinity. 

Another  problem  for  possible  future  study  is  the  case  when  the 
rate  of  the  time-optimal  control  is  bounded.  The  system  is  then  the 
following 

x»y 

y  =  -f(x)  +  z 
z  =  v  with  jv  |  <  1 , 

v  =  u  is  the  rate  of  the  control  u.  Thus  &  third  state  variable  z 
is  introduced,  and  the  new  control  v  is  bounded  in  magnitude.  The 
adjoint  system  is  then 
\  =  **'(*)  \ 

*2  *  -  \ 

=  -  Hg 

and  the  Hamiltonian  is 

H  *  T^y  +  T4,  ["*  -  f(x)^  +  v  =  h  *  constant. 

In  figure  19  *n  example  for  9  =  0.03  is  shown,  liie  integration  is 


.  > 
y 


* 


. 


' 
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performed  backward  in  time  starting  in  the  origin  with  the  following 
initial  conditions:  v  =  *-1  (the  case  v  *  -  1  is  symmetric  starting 
in  the  fourth  quadrant)  and  T)^  3  +  1«  Then  h  -  1.  Since  the  adjoint 
cyste.v  is  homogeneous  one  variable  or  tj2  can  be  chosen 
arbitrarily,  then  the  other  variable  is  determined  by  the  Hamiltonian. 
In  this  integration  the  initial  value  of  r?2  vas  varied  for  each 
trajectory.  In  figure  19  the  heavy  lines  represent  the  switching  locus 
where  v  changes  its  sign.  The  dashed  lines  mark  the  points  (x,y) 
where  z  =  0.  Much  work  has  been  done  for  linear  systems  with  bounded 
rate  of  control,  e.g.,  [7],  But  for  the  nonlinear  system  the  domain 
of  controllability  needs  to  be  determined  and  the  set  of  attainability 
should  be  described.  An  Investigation  on  the  maximum  number  of 
switchings  would  be  of  interest. 

A  further  problem  is  to  have  a  constraint  on  both  the  iragnitude 
and  the  rate  of  the  time-optimol  control.  Then  for  the  above  system 
the  control  v  and  the  new  state  variable  z  are  bounded.  Some  results 
are  available  for  linear  systems,  e.g,,  [8],  But  nothing  has  been 
done  yet  for  nonlinear  systems. 
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